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Summary. It is widely agreed upon that AGN and Quasars are driven by gas ac- 
cretion onto a supermassive black hole. The origin of the latter however still remains 
an open question. In this work we present the results of an extremely simple cosmo- 
logical model combined with an evolutionary scenario in which both the formation 
of the black hole as well as the gas accretion onto it are triggered by major mergers 
of gas-rich galaxies. Despite its very generous approximations our model reproduces 
the quasar density evolution in remarkable agreement with observations. 



1 Introduction 

While it is widely agreed upon that galaxy-galaxy interactions and, in partic- 
ular, mergers play a crucial role for the growth of supermassive black holes 
in quasars and also for providing the fuel for quasar activity, there is still 
dispute whether these black holes are of primordial nature or not. N-body 
simulations (e.g. by [1]) show that the tidal forces in interacting galaxies can 
trigger strong gas inflows towards the center of the merger. The mass of this 
gas is sufficient to build a supermassive BH of 10 7 . . . 10 10 M Q and provide 
enough fuel for quasar activity. On the basis of the /3-viscosity model by [3] 
calculations by [4] show that it takes less than 10 9 a from the merger to form 
a fully developed quasar even if no primordial supermassive BH was originally 
present. This model can be tested by comparing the resulting co-moving space 
density of quasars to those derived from observations (e.g. [5]). 



2 Outline of the model 

Based on results by [4] we assume an average time delay of 5 • 10 8 a between 
the merger and the peak quasar activity. The individual values of this delay 
depend, of course, on details like the size (sdisk) and initial mass (A/^isk) of 
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the disk. Estimates of the relevant (viscous) timescale show its dependence 
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x S disk^ r disk !• This leads to the - at first glance surprising — finding of the 
faster formation of the more massive black holes. For our present purpose, 
however, an average value of this delay serves its purpose. For more details of 
this model, we refer the reader to the contribution by Duschl and Strittmatter 
in this volume and in [4]. 

For our purpose it is sufficient to use a "test universe" comprised of 50 000 
galaxies. This test universe is expanding according to an Einstein-de Sitter- 
cosmology with a Hubble constant of Hq = 72 s ^ c . Our galaxies are treated 
as particles with a finite cross section and a thermal velocity dispersion of 
v s = 300^. In this framework we compute the merger rate for each simulated 
time bin (10 • 10 6 a) by taking into account two different processes: Direct 
geometrical hits and gravitationally driven mergers. For the first process we 
assume every galaxy to have a spherical cross section with a radius of 15 kpc. 
For the second process the cross section radius is r cs = 2G*Mv~2, with M 
being the mass of one galaxy and w re i being the relative velocity between both. 
We assume a delay between the merger and the onset of the quasar phase of 
5 • 10 8 a for the geometrical hit and 1 • 10 9 a for the gravitationally driven 
merger. The quasar phase in turn is assumed to last for another 5 • 10 8 a. 

In our model the Universe is treated in very simple manner: As an ex- 
panding spherical box. The size of this box is determined by deriving today's 
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matter density from the assumed Hubble constant ((|7rG*po) = Hq) and 
then using the Friedmann-Lemaitre-equation to calculate the according ra- 
dius at any given time. We completely neglect structure formation and start 
our simulation with large galaxies (M = 10 11 M Q ) already in place 2 • 10 8 a 
after the Big Bang. To account for the formation of galaxies we increase their 
number over the first 5 • 10 8 a of our simulation until the final number of 50 
000 is reached (see fig. 1 for the effect of this procedure). 
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Fig. 1. Co-moving space density of quasars as computed by our simulation. In panel 
a) all galaxies are in place at the start of the simulation. In panels b) and c) the 
number of galaxies is gradually increased over the first 5- 10 8 a. In panel b) a constant 
galaxy formation rate and in panel c) a Gaussian formation rate were used - the 
resulting co-moving space density of quasars is almost exactly the same. 
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3 Results 

The comparison between our model results and observational data compiled 
by [5] is shown in fig. 2. The quasar density evolution at lower redshifts and 
the position of the peak at z w 2.5 match very well. Please note that the 
X-ray selected sample contains lower luminosity AGN in addition to quasars. 
In this respect it is natural that our results resemble the 2dF curve (from 
[2]) better than the data from [5]. At higher redshifts the deviation between 
simulated and observed co-moving spatial density is increasing. In our model 
we clearly miss the earliest quasars which arise from exceptionally fast evolving 
mergers. Despite this shortcoming of our results their overall agreement with 
observational data is remarkable. 

I • ChoMP + CDF + ROSAT 2dF (Croom et at. 2004) : 

; ▲ ROSAT (Miyaji et al. 2000) □ SSG (1995) \ 
- * COMBO-17 (Wolf et al. 2003) O SDSS (Fan et al. 2001) 
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Fig. 2. Comparison of our simulated quasar density evolution with observational 
results. The broken line is the same as in panel b) of Fig. 1, while the underlying 
graph has been taken from [5]. The theoretical curve has been scaled to match the 
2dF results. 
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